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Abstract 

The utility of optical coherence tomography signal intensity for measurement of glucose concentration has been analysed in 
tissue phantom and blood samples from human subjects. The diffusion equation based calculations as well as in-vivo OCT signal 
measurements confirms the cyclic correlation of signal intensity with glucose concentration and scatterer size. 



o 

I. Introduction 

■ Everyday, almost 150 million people world wide face the problem of diabetic metabolic control. Both the hypo- and hyper- 
£X| glycemic conditions of patients have fatal consequences and warrant blood glucose monitoring at regular interval. Existing 

monitors of blood glucose can be widely classified into three classes viz., invasive, minimally invasive, and noninvasive. Invasive 
' '. monitoring require small volume of blood and they are inappropriate for continuous monitoring of blood glucose. Minimally 
O • invasive monitors analyze tissue fluid, and skin injury is minimally. On the other hand, noninvasive devices are painless and 
I skin injury is absent. 

■ Many research groups are working on to develop a real time noninvasive tool for monitoring blood glucose at clinical 
level. Recent review by Bazaev and Selishchev [1] discuss various noninvasive methods for blood glucose measurement 
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and monitoring. Existing non-invasive blood glucose monitoring work on various modalities such as absorption spectroscopy 
^ . (2), optical activity and polarimetry (3), (4), Optical Coherence Tomography (OCT) |3), 0, bioimpedance spectroscopy 0, 
. ,__| fluorescence [Tj, etc. The other optical techniques include Raman spectroscopy [8| and reflectance spectroscopy [9|. However, 
£^ due to optical interference, poor signal strength, and calibration issues, optical methods still face many challenges flOl . 
(—j The measurement of blood glucose level using OCT technique exhibits large fluctuations due to motion artifacts or other 
O ( physiological and environmental conditions. 0, @. 

1—1 Optical coherence tomography is a nondestructive technique that examines the internal structure of superficial layers of 
I biological tissues. It is based on interferometric recording of near-infrared light backscattered from the point of study, which 
£>. ■ could carry information. Conventionally, in OCT backscattered light is collected, measured, and integrated to assemble images 

o:nni. 

C*~) The impetus for the present paper comes from two reasons. Firstly, to correlate the glucose concentration with experimentally 
t^T observed diffuse reflectance signal from OCT. Secondly to examine the usability of OCT signal intensity for measurement 
. of glucose concentration. Accordingly, the paper has been divided into three parts. Section two deals with the theoretical 
analysis where Mie scattering theory and diffusion equation are used to understand the dependence of solute concentration 
on OCT signal intensity in a turbid medium. In the third section, experimental results are exhibited for various samples of 
aqueous solutions of glucose dissolved in tissue phantom as the scatterers. The experiment was also performed with blood 
• • ■ samples from voluntary blood donors. In the fourth section, the numerical calculations were carried out for various glucose 
. — concentrations in aqueous solution having different particle sizes as scatterers using the theoretical analysis. It is observed that 
a cyclic correlation of diffuse reflected signal exist as a function of the glucose concentration and size of the scatterers. To 
5— I \ support the theoretical results, experimental observations are also presented. 

II. Theoretical Formulations 

Blood is a turbid medium consisting of extracellular fluid (ECF) and various scatterers e.g. red blood cells (RBC), white blood 
cell (WBC), etc. Multiple scattering from RBC gives rise to diffusion of light from blood. The average diffuse intensity (Ud) 
is related to the scattering cross section via the relation (6) 

U d (r) = £/ d (0)exp(-fc d r). (1) 



kd is the spatial decay constant and is given by kd ( = p\lia a G tr ) where, p is the number density of scatterers in the medium, 
Otr is the transport cross section and is given by o t r [= cr s (l — g) + a a ] with cr s being the scattering coefficient, g is the 
anisotropy parameter and <j a is the absorption coefficient. 



2 



The scattering cross sections can be calculated from the Mie scattering theory as lfl2l 

2a 2 ^E„(2n+l)(-l) n ( a „-6„)| 2 



where, 



and 



bn = 



ip n (a)ip n (f3) - rmjj n {P)ij) n {a) 
rmp n (a)i>' n {fi) - ^(WnW 



(2) 



with a = k d a, k d being the propagation constant of the radiation in ECF and a is the scatterer size. /3 (= mfc^a) with m being 
the relative refractive index of the scatterer (RBC in the present case). Also, ip n {o<) = aj n (a) and £, n (a) = ah n (a), where 
j n (a) and h n (a) represent the spherical Bessel and Hankel functions, respectively. With the knowledge of these parameters 
we can calculate the path length resolved diffuse reflectance as 

1 d 



R(L S 



-U d (L)\ 



L=L, 



(3) 



U d (0) dL 

where, L s (= L s q + J2iLo -^si) is the round trip path length to the sample surface and total path length within the sample that 
accumulates during each scattering event [6|. A theoretical estimate of the OCT signal intensity can be made from the diffuse 
reflectance as follows: In an OCT set up the interference between the light scattered from the sample arm and reference arm 
is detected. The light intensity at the detector is given by 



E' s (t,L s )dL s +E r (t,r) 



(4) 



where, r[= AL/c] is the time delay corresponding to the round-trip optical path length between two beams and E' s (t,L s ) is 
the path length resolved field amplitude. E r is field amplitude from the reference mirror. 

The OCT signal obtainable from a turbid medium is a convolution of path length resolved diffuse reflectance and the low 
coherence function arising from the coherence property of the laser source. Accordingly, the OCT signal may be rewritten as, 

a/2 



I d {L s ) =I r [R{L s yi 2 ®C{L s )] 



(5) 



where, Id{L s ) is the OCT signal detected by a photodiode, I r is the signals from reference arm. C(L S ) is low-coherence 
function of the source. 

The change in the intensity of OCT signal vary in accordance with the square root of the variations in the diffuse reflectance. A 
close observation of equations (1) to (5) clearly shows that the diffuse reflectance R(L S ) is intimately related to the scattering 
cross section a s , via equation (2). Also equation (2) suggests that the scattering cross section is a function of (i) scatterer size 
and (ii) relative refractive Indies of the scatterer and the fluid in which the scatterer is dispersed. Consequently, one can expect 
a variation in the OCT signal intensity with the change in the refractive index of ECF as well as change in the size of RBC. 
The amplitude of the OCT signal A(L S ) is proportional to the diffuse reflectance (R(L S )) and is given by 

A£/ d (r) . 



A(L S 



U d {r) 



\L=L 3 



(6) 



It is worth noticing that the RBC cells, which act as source of scatterer in blood have various sizes in the subjects under study. 
Consequently, the utility of OCT setup for examining glucose concentration in blood requires the knowledge of weighted 
average of the signal (AW(L S )). Under the Gaussian distribution approximation, it is given by lfl3l 



AW{L S ) 



n £ — ' 



AU d {r) 
U d (r) 



■ exp 



S n — So 
AS 



(7) 



with ^o, S n and AS* as the size of RBC, mean size of RBC and variance in size, respectively. The reported values of RBC 
in healthy subjects are So « 7.5/im and AS « 1/im Ifl4l . 



III. Experiment 

The optical coherence tomography setup used in the present report is developed in house, using a fiber optic Michelson 
interferometer. We use superluminescent diode (SLD, Superlum SLD-371) as the light source, which is a broad band source 
with a center wavelength of 841 nm and line width of 52nm. The selection of wavelength in the near infrared region has the 
advantages of large scattering coefficient and low absorption coefficient in tissues. However, these coefficients are altered when 
index mismatch occurs between ECF and cells. 
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Figure 1. Schematic of the Optical Coherence Tomography setup. SLD - Superluminescent diode, DC - 3dB bi-Directional Coupler, PD-Photo Diode, 
S - Sample, R-Reference arm with scanning assembly, BPF-Band pass filter, ED-Envelope detector. 



A. Sample Preparation 

The primary goal of the present work is to find the correlation between the amplitude of OCT signal with concentration of 
glucose. A systematic measurements were carried out with the following samples; (i) tissue phantom based as intralipid lfT31 
and (ii) blood samples from voluntary human subjects. 



B. Tissue Phantom 

Intralipid is an emulsion of soy bean oil, egg phospholipids and glycerin [16|. Intralipid is widely used in optical experiments 
to simulate the scattering properties of biological tissues and can be use as a good scatterer as RBC in blood. The major 
advantages of intralipid are its well known optical properties and the similarity of its microparticles to lipid cell membranes 
and organelles that constitute the source of scattering in biological tissue ifTTIl . We use intralipid as a tissue phantom that 
provides the backscattered component. Average size of scatterers in intralipid measured using confocal microscope is found to 
be 3.5 /im. 

Using standard approach, we generate glucose concentrations of - 200 mg/dl in 1ml sample tissue phantom for experimental 
observations. This range covers hypoglycemic (< 80 mg/dl), normal (80-120 mg/dl) as well as hyperglycemic (>120 mg/dl) 
conditions. A settlement time of 2 minutes was given after addition of glucose to the intralipid solution so that the added 
glucose could alter the scattering properties of the medium. Simultaneously, the change in refractive index of the solution is 
also monitored using Abbe's refractometer. Measurement was done using fixed volume 1ml of tissue phantom and diluted to 
0.01% in distilled water (100ml) and concentrated glucose solution prepared by dissolving 200 mg glucose in 100 ml distilled 
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water. The concentration of glucose in intralipid is increased in steps of lOmg/dl using U-TEK Chromatography syringe with 
least count of 5/il. 



C. Blood samples 

Blood samples are collected from Pathological Laboratory. In order to optimize the setup with human subjects, blood samples 
were collected from various voluntary subjects. Conventional pathological method known as GOD/POD (Glucose Oxydised / 
Peroxydised) using semi-automated blood analyser was used simultaneously to measure the exact value of blood glucose. A 
smaller part (fewfil) of the glucose solution was used in the OCT setup. Small quantity of 4 /A blood taken using a Eppendorf 
Pipette. 



IV. Results 

To confirm the changes in the optical properties of the solution with increasing glucose concentration in intralipid, Abbe 
refractometer was used to obtain the refractive index of tissue phantom and found to be 1.407. As discussed earlier, addition of 
50/iZ glucose to intralipid solution of volume 1ml increases the glucose concentration by lOmg/dl. At the same time, it reduces 
the index matching between the scatterer and the medium. The measurements of refractive index and scattering coefficients 
using Abbe refractometer and OCT, respectively are made simultaneously to avoid the temporal fluctuations in values. Results 
obtained with refractometer are displayed in Fig. [2] A linear relation is obtained between the glucose concentration and the 
refractive index. The data from Fig. [2] is used to obtain an empirical relation between the glucose concentration and the 
refractive index. 
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Figure 2. Variation of refractive index of intralipid with Glucose concentration measured using Abbe's refractometer. 



OCT signal amplitudes obtained with tissue phantom and for various glucose concentrations (10-200mg/dl) are shown in Fig. [3] 
In order to reduce the noise in OCT signal, we displayed the normalised average area of scattered signal within the coherence 
length and shown as dots. The OCT signal amplitude as obtained now is numerically equal to the average diffuse intensity Id 
defined in eq. The theoretical fit shown as solid line in the Fig. [3] is the weighted average as discussed earlier. The results 
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indicate an oscillatory nature of backscattered component of the OCT signal with glucose concentration. These observations 
with tissue phantom show a cyclic correlation between glucose concentration and OCT signal amplitude and are in good 
agreement with the theoretically calculated values of diffuse reflectance (solid line). Both the curves exhibit similar nature 
which confirms the usability of present theoretical calculations to correlate them with blood glucose levels in human subjects. 
In order to confirm the same with clinical measurements with human subjects, we collected samples of voluntary donors from 
pathological laboratory and the results are discussed below. 

Simultaneous measurements of blood glucose was carried out using chemical methods as described in previous sections. The 
correlation between blood glucose level of human subjects and the amplitude of the OCT signal are represented in Fig. [4] The 
curve definitely shows a cyclic correlation, similar to the results predicted from eq. (|4]i and Fig. [3] 




Glucose Concentration (mg/dL) 

Figure 3. OCT signal amplitude obtained in a sample of tissue phantom with different glucose concentrations are shown as filled circles with error bars. 
The solid line is the weighted average obtained theoretically using equation {7}. 

We therefore, find that the cyclic correlation of glucose concentration with OCT signal restricts the OCT signal intensity based 
glucose measurements in human subjects. 

To conclude, although it is well known that the optical properties of tissue vary with glucose concentration and the index 
mismatch between the ECF and scatterers. There are many factors which will affect the OCT signal intensity. For example 
the refractive index of ECF will change with change in any types of soluble impurities (example : urea) and the change in 
the scatterer size. Accordingly, well defined predictions for glucose concentration measurements using OCT signal intensity 
are not possible because of the non-uniform size distribution of various particles as well as presence of creatinine, lactic acid, 
serum albumin, and NaCl, Urea, etc. in blood. Therefore, the results presented here may not be used as a tool for estimation 
of blood glucose levels, however, it throws light on the understanding of light scattering from blood samples, in the presence 
of various types of scatterers differing in physical and optical properties. 
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Figure 4. The filled circles represent the amplitude of OCT signal obtained from various blood samples of voluntary donors having different glycemic levels 
while the solid curve is the weighted average obtained theoretically using equation 0. 
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